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Abstract. We consider a two dimensional semiconductor with a local attraction among the carriers. We 
study the ground state of this system as a function of the semiconductor gap. We find a direct transition 
from a superconducting to an insulating phase for no doping at a critical value, the single particle excitations 
being always gapped. For finite doping we find a smooth crossover. We calculate the critical temperature 
due to both the particle excitations and the Berezinkii-Kosterlitz-Thouless transition. 



PACS. 74.20.-z Theories and models of superconducting state - 
- 74.20.Mn Nonconventional mechanisms 



74.20.Fg BCS theory and its development 



1 Introduction 



In the standard BCS theory of superconductors the gap A 
in the spectrum of quasiparticles is due to the breakdown 
of gauge symmetry at the phase transition, in much the 
same way as the gap in a spin density wave state is due 
to Bragg scattering off the doubled periodicity. The gap 
disappears with order and consequently A — 1.76 Tc is 
comparable to the transition temperature Tc. Largely mo- 
tivated by high temperature superconductivity, interest 
has grown recently in a "pseudogap" situation, in which 
a somewhat blurred gap A persists in the normal state, 
such that A ^ Tc Q. There are many models for such a 
pseudogap, the simplest one being a transition from the 
BCS limit towards a "Bose Einstein condensation of pre- 
formed pairs" : if the attraction between electrons is strong 
enough, bound pairs (let us say singlet) form first, with a 
binding energy eb, then below a critical temperature Tc 
they Bose condense. It was shown long ago that the evolu- 
tion between the two limits is smooth, with however a very 
different physics in the strong binding limit M . Then the 
gap A is just the binding energy, a "molecular" quantity 
which has nothing to do with symmetry breaking; in con- 
trast Tc is controlled by center of mass motion of bound 
pairs, which are not broken: as a result A ^ Tc, a simple 
minded interpretation of the pseudogap. A similar evolu- 
tion occurs for the transition from a paramagnetic metal 
to an antiferromagnetic insulator in half filled band: for a 
strong repulsion the Mott gap A exists irrespective of the 
magnetic order. The latter is only a secondary feature, 
while it was a key actor in the formation of A at weak 
coupling M. 
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In order to build an accurate description of realistic 
systems one needs an interpolation between these two lim- 
its. That means including quantum fluctuations of the or- 
der parameter A, (mostly of its phase) which are system- 
atically ignored in a mean field BCS calculation. That is 
a difficult task and the literature on the subject is large 
[||. We will not attempt to review it: as usual in perturba- 
tion theories it relies on approximations that are not re- 
ally controlled. We will just make a few simple qualitative 
comments. Grossly speaking one expects the crossover be- 
tween the two limits to occur when the (zero temperature) 
gap A is comparable to the Fermi energy ep'- they are the 
only energies in the problem. That can be checked explic- 
itly for a two dimensional system, in which phase fluctu- 
ations are enhanced, leading to the well known Kosterlitz 
Thouless transition [g| . The resulting transition tempera- 
ture T^^ is controlled by the phase stiffness of the ground 
state: it can be compared to the mean field Tf^^ which 
accounts only for pair breaking. Such an elementary calcu- 
lation is sketched in section |j: it shows that the crossover 
T^^ ^ j-BCS occurs when A « £p jh. There a bound pair 
would sit in the middle of a continuum, which makes no 
sense. Sure in the very strong coupling limit one returns to 
preformed bound pairs, but there exists a broad interme- 
diate region near the crossover where phase fluctuations 
are already dominant without any bound state. 

In the above model the same interaction is responsible 
for the gap A and for Bose condensation. That may well 
be the case and the future will tell whether high tempera- 
ture superconductors pertain to such a mechanism or not. 
However that is not a necessity: one may also imagine that 
a primary gap Ao is due to some other mechanism that 
has nothing to do with superconductivity. That may be 
for instance a lattice distortion of some sort or a magnetic 
instability, as suggested by various authors g : the start- 
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ing point is then a semiconductor rather than a metal. If a 
BCS attraction between electrons is added to such a sys- 
tem, it may lead to superconductivity, even if there were 
no free carriers to start from: free pairs appear sponta- 
neously if the gain in binding energy is larger than the cost 
to produce free electrons and holes across the gap. Such a 
mechanism is actually familiar: it was introduced long ago 
by W. Kohn in his theory of excitonic insulators . There 
the Coulomb interaction is repulsive and the bound pair 
(the exciton) involves an electron and a hole: excitons ap- 
pear spontaneously if their binding energy is larger than 
the gap. What we present here is just the reverse situation 
in which electrons attract: bound pairs involve two elec- 
trons or two holes; at T = they Bose condense (bound 
carriers are then time reversed of each other) , hence super- 
conductivity. It is clear that superconductivity will modify 
the real gap Ag - hence the name "primary gap" for Aq. 

The possibility of superconductivity in a semiconduc- 
tor is not new: it appeared since the beginning of the BCS 
"era" , as emphasized by the review article of M.L. Cohen 
1^ (for subsequent development, see ||]). However only 
intraband interactions in a doped material were included, 
the only new feature being the scarcity of carriers. Here we 
are concerned with the effect of interhand interactions, in 
which carriers in both bands conspire in establishing su- 
perconducting order. That was also considered before |10| , 
but in a somewhat different context of interband pairing, 
while we consider only the usual intraband pairing of time- 
reversed states. We show in Appendix H that the latter is 
energetically preferred. (An obvious question is of course 
whether they can occur together.) 

Many features of such a "superconductivity in a semi- 
conductor" are general, in the sense that they do not 
depend on the physical mechanism responsible for gap 
opening. We fell that such an interplay of Ao and su- 
perconducting order deserves study, the more as it leads 
to somewhat unexpected simple conclusions. The present 
paper offers a very naive discussion of that problem. At 
this stage we have no pretense whatsoever to explain any 
realistic material: our only purpose is to identify and to ex- 
plain qualitative features of the various limits at hand. In 
order to be enlightening and hopefully convincing, such an 
analysis must be simple: we therefore choose deliberately 
the most oversimplified model and formulation, in such a 
way as to have transparent results. We are perfectly aware 
of the many loopholes and traps of the model, and we will 
try to point those which we feel are important. Starting a 
discussion of these many complications would however be 
meaningless at a stage where even the feasibility of such 
a mechanism is an open issue. 

In section we introduce our model, and we discuss the 
simplest situation in which the semiconductor is undoped: 
were it not for superconductivity it would be an insula- 
tor at T = 0. Within a mean field BCS approximation 
the algebra is elementary and it leads to a direct transi- 
tion from insulator to superconductor at zero temperature 
(when the primary gap Aq reaches a critical A*^). At finite 
T one easily calculates the order parameter Z\, the actual 
gap Ag and the critical temperature Tc : as could be sur- 



mised a pseudogap situation prevails near Z\*. Section rais 
concerned with slightly doped systems: then the ground 
state is always conducting. With an attraction between 
carriers it is always superconducting. Still within a BCS 
mean field approximation, we explore the various limits. 
The interesting case is Ag > A*, where everything is due 
to doping; the behaviour of gaps and critical temperature 
in that limit is somewhat unexpected. Up to that level 
the dimension of space and the nature of the gap in the 
Brillouin zone do not matter much: they become crucial 
in section which is concerned with a much more deli- 
cate issue, namely phase stiffness and Kosterlitz Thouless 
fluctuations. The analysis is here more speculative: as far 
as we can see the transition is always monitored by phase 
fluctuations in realistic models for the gap, in contrast to 
ordinary BCS metals. The reason is a very small phase 
stiffness, which would vanish for a filled band. In a brief 
conclusion we emphasize the many questions that remain 
open. 



2 The BCS model for an undoped 
semiconductor 

We start from a normal Fermi liquid that has a gap 2Z\o, 
and we choose the origin of energy at midgap. The gap 
may be direct or indirect in momentum space: at that 
stage it does not matter. Let p{^) be the density of states 
(per spin): in a first approximation we take p{£^) con- 
stant: everything is symmetric around ^ = 0. Since the 
total number of states must remain the same, opening 
the gap Ao means that the conduction band shifts from 
[0,Wm] to [Z\o,^o + Wm] while the valence band goes to 
[—Ao — LOrm —Ao]- The conduction and valence bands are 
just shifted away from each other. tOm is a bandwidth 
which we assume ^ Ao (the gap is small). The choice 
of a constant p is in accordance with our choice of sim- 
plicity, but it is not realistic for two reasons (i) at band 
edges p goes as ^^/^ in three dimensions, {ii) in usual sit- 
uations states that are repelled off the gap remain in its 
vicinity, in an energy range which has the same order of 
magnitude Ao (see Fig. |]); here these states are repelled 
at infinity. We will show in Appendix |A| that more realis- 
tic choices for p{(_) do not change the results qualitatively: 
they are not worth the price of complications. 

We then add a local attraction U between fermions 
and we look for a BCS superconducting state in which 
electrons (or holes) are paired in singlet states with zero 
total momentum, (k t) with (— k |) (one the time reversed 
of the other). If the bottom of the band is made up of 
several valleys at band edge, that means pairing of elec- 
trons between symmetric valleys. We treat the problem 
within a mean field approximation: the Hamiltonian con- 
tains anomalous terms 



A CkT c_k 



ti 



where A is the usual order parameter, not to be confused 
with gaps, which must be determined self-consistently. Be- 
cause of electron- hole symmetry the chemical potential p 
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Fig. 1. Possible modification of the density of states near the 
band edge. The thick line correspond to the approximation 
used in the text. The dotted line is instead an example of a 
more realistic 2d case where the removed stated remain near 
the edge. The dashed line is a typical 3d case where the density 
of states vanishes at the edge with a square root dependence 
on the energy. 
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Fig. 2. Undoped case, we report as a function of Aoj Am the 
following quantities: the order parameter Aj Am (continuous 
line), the excitation gap AgjAm (dashed line), the critical tem- 
perature Tel Am due to pair breaking (long dashed line), and 
the Kosterlitz-Thouless critical temperature T^'^/Z!™,, [see 
section B (dotted line). 



of our undoped system remains at midgap and the reduced 
kinetic energy ^^ = [e^ — A*) starts at -^Ao instead of 0. 
The algebra is identical to that of BCS except for that 
shift in the energy integration. In the absence of a gap, 
Ao = 0, the order parameter A would take the regular 
metallic value Am given by the famihar BCS self consis- 
tency condition 



pU_ 
2 



-l-l^r. 



di 



^e + ^. 



pU log 



2w„ 



It is convenient to characterize the attraction by Am in- 
stead of U. Note that for a regular metal the supercon- 
ducting gap is just Ag = A. In the presence of a semicon- 
ductor gap Ao the self consistency equation for the new 
order parameter A becomes 



1 = pU log 



2uj„ 



Ao + y/Al + Z\2 



(2) 



In order for the two equations to have a solution with the 
same U, we must have Am = Ao + \J ^o + ^ i 



I.e. 



A= ^jAm {Am - 2Z\. 



(3) 



Superconductivity persists with a reduced order parame- 
ter until the gap 2Z\o reaches the metallic BCS gap Am'- 
then the gain in interaction energy due to superconduc- 
tivity can no longer overcome the cost of kinetic energy in 
producing free carriers across the gap. Note the complete 
analogy with excitonic insulators in which condensation 
produces electron-hole pairs in the presence of a repulsive 
U 0. Such a model provides an elementary example of 
a direct superconductor-insulator transition, which here 



is sharp at T = 0. The quasiparticle gap Ag is just Ao 
in the insulating state {A = 0). In the superconducting 
state it is \/Al + A? , i.e. Ag = Am — Ao. The situation 
is illustrated on Fig. 0: note that the gap is finite when 
superconductivity disappears. 

The calculation is easily extended at finite tempera- 
tures. The order parameter A becomes a function of T. 
We need only introduce the fermionic occupation factors 
in the gap equation: 



d^ 



tanh/3Vr+^/2 

VeTA^ 



d^ 



ve+^. 



(4) 



(we retain the zero temperature Am as a measure of the 
attraction U). The integration path P takes into account 
the absence of states in the middle of the band: P = 
[-{Ao+LOm),-Ao] and [+Ao,+{uJm + Ao)]. For Z\o = 0, 
Eq. (Q) is the standard BCS gap equation, yielding the 
usual A{T) and Tc. At the other end the critical temper- 
ature goes to zero at the superconductor-insulator transi- 
tion Z\* = Z\m/2 (there Z\ = at zero temperature: there 
is no need of a T^ to make it vanish). The behaviour of 
Tc{Ao) is shown on Fig. 2. 

The approach to the transition is interesting: we sub- 
tract the Ao < A* Eq. (0) from the same equation at the 
transition, Ao = A*: 



+ 00 



1- tanh/3cC/2 



d^- 



A* 



tanh/3cC/2 



d^ = 0. (5) 



When Ao -^ A*, Tc vanishes and we can approximate 
tanha; with 1 — 2e^^^. The first integral reduces then to 
the Exponential integral function and for large values of 
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Fig. 3. A{T) for r = 1, 0.9, 10"\ 10"^ and 10"^ The evolu- 
tion from the familiar BCS shape for large t to a new one for 
r small is evident. 



(3cAl we obtain: 



Tr_ 



Al A 



Al-A, _ 

— T . 



(6) 



Defining the function w{x) as the solution of the tran- 
scendental equation x = we^ , we can write Tc in terms 
of w{x) 

A* 

w(2/tj 

The function w{x) is a generalization of the logarithm and 
for x 3> 1 has the following expansion: 



'w{x) = In .T — ln[ln a: — ln[ln x — In . . .]] , 



(8) 



The convergence of that expansion is very slow [shown 
in (^], but the first term is enough to provide a rough 
result: the critical temperature vanishes logarithmically 
inT[T,~Z\:/ln(2/r)]. 

Altogether the presence of a semiconductor gap Ao 
acts to decrease T^, a result which is reasonable. But such 
a decrease should be compared to that of the zero tem- 
perature order parameter Zi, which behaves as ^/t: we 
see that the ratio A{T = 0)/rc becomes extremely small 
when Ao -^ A*^, in contrast to the usual BCS result. The 
physical origin of that behaviour is clear: due to the finite 
gap Ag the quasi particle excitations are much less effi- 
cient in destroying the superconducting order than they 
would be in the BCS case, and consequently Tc increases. 
Such a behaviour is exemplified in Fig. H in which we plot 
the order parameter A{T) [obtained from solving (|4|)] for 
various values of Ao/Am'- the elongated low plateau near 
the transition is striking. 

All that discussion is based on a BCS like mean field 
approximation, in which condensed pair breaking is the 
only allowed thermal excitation, ultimately responsible 
for the destruction of superconductivity at T,.. Such an 
approximation is by no means obvious and it should be 



kept in mind. In practice there exist other excitations that 
might be more efficient in destroying superconductivity. In 
particular, phase fluctuations are gapless excitations and 
at least for Aq near Z\* they are expected to reduce dras- 
tically the critical temperature. We will discuss this point 
more extensively in section for the two-dimensional sys- 
tem. There we find that phase fluctuations seem to be 
always responsible for the transition to the normal phase, 
even for Ao = 0, a somewhat unexpected conclusion. 



3 Doped system 

We now extend our problem, allowing for a flnite dop- 
ing. We do it because it leads to new interesting physics: 
the sharp transition at zi* becomes a smooth crossover, 
and superconductivity persists at zero temperature for any 
value of Ao- Beyond that crossover, the superconducting 
behaviour is very unconventional over a broad range of pa- 
rameters, in contrast to the undoped case in which queer 
behaviour occurred only at a specific transition point. In 
that respect the doped system is much more generic. 

Starting from the non interacting, metallic system at 
T = we suppose that the upper band is populated up to 
an energy /i„j = e, referred to midgap. When the gap Ao 
is restored, the "semiconductor" is doped by an amount 
5N = 2pe (referred to the filled valence band). In the 
metallic state Ao — 0, chemical potential and supercon- 
ductivity are unaffected by each other (the BCS picture is 
symmetric with respect to Fermi level): /j,™ remains equal 
to £. In the opposite limit of an insulator with no supercon- 
ductivity, /i — Ao+s is pushed into the conduction band, a 
familiar feature of degenerate semiconductors. In between 
we expect no sharp superconductor-insulator transition: a 
weak superconducting order due to the few free carriers 
in the conduction band should persist at any Aq. Keep- 
ing Am and e fixed, we thus expect the zero temperature 
chemical potential fi and order parameter A to depend 
on Ao as showed respectively on Fig. |j (a) and (b). The 
transition around Ao = Am/'^ is a crossover, but with a 
large change of /i. 

Let us assume an arbitrary temperature T. In order 
to fix the two parameters /i(T) and A{T) we need two 
conditions. One is the self consistency equation for the 
order parameter, which is the same as (0) but for the 
introduction of the chemical potential in the expression 
for the quasiparticle energy: 



tanh/3£:(0/2 



dc- 



rfe 



ve+^ 



(9) 



where E{S^) — ^/(C ^ m)^ + A"^. The other one is the con- 
servation of particle number, stating that SN is unchanged: 



2p I d( 



vHo + fF[Em 






2p 



(10) 
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Fig. 4. (a) Chemical potential (/i) for different value of ej Am 
as a function of Aoj Am- The dashed line is the bottom of the 
upper band (^AojAm)- (b) The same for the order parameter 

AjAm- 



3.1 Ground state properties 

We begin by considering the zero temperature limit, 
this case integrations are elementary and we obtain: 



In 



Z\2 = 



A„ 



(zio - /i)' + A^ 



Ao^[i 



(Z\o + /i)' + /i2 



_, (12) 

System ( p^ ) fixes /i and Z\ as a function of the gap Ao 
and doping bN = 2p£. Once they are known, the overall 
gap for elementary excitations is either A {\i fi > Ac) or 



{Ao — /i) + A"^ (if ^ < Ao). Numerical solution of these 
equations is easy and it is shown in Fig. |4[ Simple limits 
can be assessed analytically as long as £ <C Am, shedding 
some light on the physical meaning of these results. 

Consider first the "conducting" side. Here free carriers 

extend deep into the valence and conduction bands: the 

chemical potential for small doping will remain close to 

midgap. We can therefore expand the second equation of 

HI) into 



£ = ^ 



1- 



A„ 



^Al + Z\2 



(13) 



Doping is a small perturbation and we can use the un- 
doped result for A: wc thus find 



^111. ^o 



Ai = 



Ar. 



2Ao 



£ + 0(£2). 



(14) 



The approximation breaks down when £/Z\* '^ t. The 
crossover is very close to the superconductor-insulator sharp 
transition found when 8N = 0. Substituting that first or- 
der result for /i into the first of (O) we find the correction 

to A: 



where /i.(^) = l/(e''^ - 
temperature /3 and 



1) is the Fermi factor for inverse 



v\i)^l~u\i) = \ 



e-M 



E{i) 



(11) 



is the familiar BCS occupation factor. The cut-off LOm can 
be sent to infinity in Eqs. (P)-(p^): we are left with a sys- 
tem of two coupled equations. 



^ = 2^:v^+(^) ^/i: + o(£4). (15) 

As long as Ao < A* doping is a minor perturbation. 

In the opposite limit, we first note that an insulating 
state at T = is impossible. It would mean a, A = 
solution to the gap equation, i.e. Af^ = 4(Z\q — /x^), which 
in turn implies /i < Aq. But with no A that means 6N = 
0. Thus superconductivity must either extend all the way, 
or display a first order transition: a sharp second order 
transition to the insulating state is precluded. We now 
show that the present model has a smooth crossover. 

Deep in the insulating limit, Aq > A*, the role of equa- 
tions is interchanged. The order parameter A will turn out 
to be small: consequently the gap equation which used to 
fix A now fixes /i 



<Al^^,') 



^j- 



Al 



(16) 
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Such an expansion is valid as long as {Ao — /i) 3> /i, which 
implies 

^^»A (17) 

fj, is still inside the band gap Aq. Granted ^, the other 
equation for SN fixes the order parameter A 



SN 



2A^ 

A '' 



Al-A* 



^ A = 



A,. 



[Al - Af] 



21 1/4 



(18) 

As surmised A is small for small doping e, which justifies 
our calculation a posteriori. The condition (n^ now reads 



£< 



8zio 



(19) 



Doping should remain small (or the gap Aq not exceed- 
ingly large). 

When the gap exceeds that range the expansion in 
powers of A is no longer warranted. But as long as doping 
£ is small compared to Z\* we expect the chemical poten- 
tial to remain close to Z\o, either below (inside the gap) or 
above (inside the continuum). We may then assume safely 
that {Ao + /-i) ~ 2Ao » A. The two basic equations ( [12|) 
simplify into 




Ao-^l + ^{Ao-^ir + A^ 

Il-Ao + ^{Ao-iif + A^ ' 



Ao 



(20) 



Multiplying one by the other we obtain the very simple 

result 

2j^/\*2 
A^ ^ ° 
Ao 

which joins smoothly with our former expansion dlq), and 
which remains valid throughout the crossover \Ao — /i| ~ 
A. Similarly we may eliminate the square roots from these 
two equations, which yields the chemical potential /i 



^1 = Ao 



A 



*2 



- e — 



2Ao 



(21) 



For very large gaps Ao ~> A*^ je, the last term is negli- 
gible and we recover the usual result for a doped normal 
semiconductor, /j, = Ao \- e. 

Let us explore these two regimes in some more detail. 
When the chemical potential lies inside the gap Ao the 
real gap Ag is just [Ao — fi) ^ A, of order Z\*^/2Z\o as fi 
approaches the band edge. The relevant ratio is 



A 

a: 



ISsAo 



A* 



<1 



(22) 



Note the strong similarity between such a regime and the 
Bose condensation of preformed pairs |2|. In the latter 
case the chemical potential fi is close to the energy of 
bound pairs, £6/2 per fermion and therefore it lies inside 
the gap. Here the reduced kinetic energy (^ — fi) likewise 



never vanishes: (Ao — /i) plays the role of efc/2. The result- 
ing behaviour is similar in both cases, with a completely 
different physics. Here the hierarchy of energy scales is 

e < Z\ < Z\g 

(indeed one has Z\^ = AeAg). In the opposite limit the 
chemical potential enters the conduction band: then the 
real gap Ag is just the order parameter A. The relevant 
ratio is that of A to the Fermi energy e measured from 
the bottom of the band 



A 

e 



l2Af 
eAo 



Except for a numerical factor 4 we see that /I <C e: we 
recover an ordinary BCS regime in which superconducting 
pairing is a small perturbation. 

What is unusual is the behaviour of the order param- 
eter A K, -^e in that "insulating regime" , for small doping 
e and moderate semiconductor gap Ao- Such an unusual 
dependence A « y^ is due to the interplay of doping and 
superconductivity: (i) extrinsic carriers are needed in or- 
der to start superconductivity, but {ii) the latter does add 
further carriers, thereby modifying A. That result for the 
dependence on doping away from Zi* should be compared 
with the result A = 2Z\*-yT for the undoped system close 
to Z\*. The two variables r and e act in the same way in 
driving the system away from the T = critical point (we 
will find this behaviour again for the critical temperature). 

Analytic calculations are possible only in limiting cases. 
In the crossover region one must resort to numerical solu- 
tions, which are quite straightforward. The result is shown 
in Fig. y where system (y^ is solved numerically for the 
values of e indicated. Both A and /x clearly display the 
transition from a "metallic" to a "semi insulating" regime. 
In the latter, /i crosses the band edge when e ~ A*^ /2Ao, 
separating Bose Einstein like and BCS like behaviours. 
Note that the chemical potential crosses the band edge 
Ao twice when e is smaller than a critical value (see also 
the discussion after Fig. ItI). 



3.2 Finite temperature behaviour 

Having a complete view of the zero temperature doped 
system, we consider now the effect of a finite temperature. 
In particular we want to investigate the critical temper- 
ature of the system. Setting A = Q into Eqs. (|9|)-(|l0|) we 
obtain: 
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Fig. 5. Critical temperature Tc as a function of Ao/Am for 
e/A^ = 0, 10-^ 10-^ 10-\ and 1. 



These equations can be easily solved numerically and the 
results are shown in Fig. for different values of e includ- 
ing the undoped case. As expected from the zero tem- 
perature analysis, the critical temperature also exhibits 
a smooth crossover from a large to a small value, and it 
never vanishes for finite doping. The region Aq < A* is 
not much different from a standard BCS superconductor: 
a comparison of Fig. || with Fig. || (a) clearly shows that 
Tc is proportional to A{T = 0) (with the same BCS factor 
for A^O). 

A new and interesting behaviour is instead found when 
Ao > zi*. In this region the critical temperature appears 
to be much larger than the order parameter A{T — 0) 
when e/ Am <C 1 . This point can easily be verified ana- 
lytically. In this limit we can substitute the ground state 
chemical potential ^o in Eq. (p4 ). Since the critical tem- 
perature vanishes for e = while Ao — /Xq > the ex- 
ponentials in Eq. (124) are large and we can expand the 
expressions. At leading order in exp{(/io — Ao)/2Tc} we 
obtain the following equation for Tc'. 



Tc^ee^^^. (25) 

Again the solution can be written in terms of w{x): 

T, = ^iiZli^ ^ 4^iZif^ . (26) 



e 


In 


e 



It is then clear that the ratio Tc/A{T = 0) diverges for e 
going to zero. One should note that the behaviour found 
for both A{T — 0) and Tc as a function of r is exactly the 
same of that found for these two quantities as a function 
of e. The origin is the same since the gap at Tc in this case 
becomes Ao — n- 

We report in Fig. o the ratio between the critical tem- 
perature Tc and the other two relevant energy scales, the 
zero temperature order parameter A(T = 0) and excita- 
tion gap Ag{T — 0) for different values of e/Am going 
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Fig. 6. Tc/Ag{T = 0) (continuous line) and Tc/A{T = 0) 
(dashed line) for different values of e/Am (indicated on the 
plot). The lines corresponding to e/Am ~ 10~^ are nearly con- 
stant at the BCS value of ~ 0.57. 



from to 10^^. For small Ag/Am the familiar BCS re- 
sults of e'' /tt fa 0.57 is recovered, while for Ao > 4im/2 a 
totally different behaviour is found. The critical temper- 
ature appears to be in between the other two quantities 
(for e < Z\„0 



Z\(T = 0) <Tc< Ag{T = 0). 



(27) 



This can be easily checked analytically, since for s going 
to zero Ag stays finite [Ag -^ Aq — \/Al — Z\^/4], while 
Tc vanishes logarithmically and A{T — 0) vanishes like 
y/e. Thus at mean field level we have found a rather pecu- 
liar form of superconductivity, whose behaviour is totally 
different from the BCS one. We note that the BCS be- 
havior is recovered also for large values of Ao when the 
chemical potential lies in the conduction band. In Fig. 
all the lines will join the BCS value 0.57 for fixed e and 
large enough Ao- If £/^m is small enough, the parameter 
Ao drives the system through four different regimes. For 
very small Ao the system is in a two-band BCS regime. 
Increasing Ao the chemical potential enters the gap and 
the smoothed transition from a normal superconductor 
to a "semi-conductor" superconductor takes place. When 
Ao ^ zi* and the chemical potential is still in the gap, 
the carriers in the upper band act as a strong coupling 
superconductor, while the lower band is nearly decoupled. 
For even larger Zio(> ^o^/2e) the chemical potential lies 
again in the valence band and a BCS behaviour is recov- 
ered for the upper band. The schematic phase diagram is 
plotted in Fig. [7| where the region of positive and nega- 
tive \i — Ao are identified on the plane Ao-e. The path 
described above correspond to the dashed line. 

All that discussion is based on a mean field, BCS like 
formulation. We should now ask whether such an approxi- 
mation is warranted. In the next section we inquire about 
the effect of phase fluctuations in two dimensions, which 
are known to provide an alternate mechanism for the crit- 
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Fig. 7. Phase diagram for the chemical potential. The shaded 
part is the region with /i < Ao- The thick line represents the 
insulating region of the phase diagram, for any other value 
of £ and Ao the system is superconducting. The dashed line 
represents the path in the phase diagram described in the text. 



ical temperature Tc, based on the creation of vortices in- 
stead of pair breaking. 



4 Phase stiffness and Kosterlitz— Thouless 
fluctuations 

We limit our discussion to a two dimensional system for 
which they are known to be crucial, leading to the so called 
Kosterlitz-Thouless transition ||^. Although there is no 
strict long range order in 2d, the low temperature decay 
of the order parameter is algebraic up to a sharp critical 
temperature T^'^ at which the decay becomes exponen- 
tial: such a transition from a "quasi superfluid" phase to 
a normal phase is due to the unbinding of vortex pairs. 
In the most naive approximation one may assume that 
the modulus of the order parameter is constant, the only 
relevant variable being its phase '?!>(x). The corresponding 
fluctuations are controlled by the phase stiffness J defined 
by the following term in the Hamiltonian 



AH 



J 



(fx\V(l){x)\'^ 



(28) 



J measures the resistance of the system to phase gradi- 
ents. From J one infers in 2d the well known exact critical 
temperature [0 



T. 



KT 



-l^' 



(29) 



The stiffness that enters (^9|) is actually temperature de- 
pendent, J{T). That dependence occurs both at the mean 
field level, through pair breaking thermal fluctuations, and 
via thermal fluctuations of the phase itself, which can be 
described within a renormalization scheme. We disregard 
the latter and we discuss briefly the mean field J^''^{T) 
which provides a neat picture of the competition between 



BCS and phase fluctuations [0|. J*^^ vanishes at the 
mean field transition temperature T^^^, hence two pos- 
sibihties: (i) J^^^(O) » T^. phase fluctuations are essen- 
tially negligible except very close to Tc- then T^'^{T) be- 
comes comparable to T and BCS turns into a Kosterlitz- 
Thouless regime; the ultimate transition is deflnitely in 
the Kosterlitz-Thouless universality class, but the rele- 
vant range is narrow and the real Tc is very close to 
T*^^. Basically BCS is correct, [n) J*^^(0) < T^. then 
the Kosterlitz-Thouless transition occurs well before pair 
breaking fluctuations become relevant: but for renormal- 
ization corrections Tc is T^^{Q). 

What we need is an estimate of the stiffness J(0). That 
is easily achieved considering a superfluid flow in which 
pairs are condensed with center of mass momentum Q/2 
instead of 0. The corresponding phase of the order pa- 
rameter is = Q • r and its gradient is Q. In order to 
calculate the energy of that flow we need a specific model 
of the band structure in momentum space: the knowledge 
of the density of states is not enough. Let us first consider 
the standard model of a 2d free electron gas, with no gap 
at all. The Fermi surface is a sphere with radius kp^ the 
density of electrons with either spin is iV = kp/2TT. For 
a parabolic dispersion law, with mass m, the center of 
mass and relative motion of a pair are completely decou- 
pled: the cost in kinetic energy of the superfluid flow is N 
fi^Q^/Sm. The phase stiffness is consequently 



J = 



Am 



477 



The corresponding transition temperature is T^'^ — ep/8, 
comparable to the Fermi energy, but noticeably smaller 
due to the accumulation of numerical factors. That must 
be compared to the BCS T^ w .57 A; as mentioned in 
the introduction the crossover occurs when A is about a 
fifth of £j?, making the idea of a preformed bound state 
meaningless. Usual low temperature superconductors are 
way below that crossover and phase fluctuations are ac- 
cordingly negligible. 

We now turn to our gapped semiconductor: we need to 
specify the detailed form of the two bands. The simplest 
guess assumes valence and conduction bands that have a 
well defined extremum in the Brillouin zone, whether at 
the same place (direct gap), or at different places (indi- 
rect gap). Near the band edge the dispersion is parabohc, 
characterized by two (different) masses (the values of these 
masses will not enter in the final expression of the criti- 
cal temperature) . In a 2d system the density of states at 
band edge is finite, and the assumption of a constant p 
is not crucial: up to that point, our model raises no diffi- 
culty. It does however when we try to calculate the phase 
stiffness J. Before attempting explicit calculations, let us 
explain the point qualitatively. In order to set up a su- 
perfluid flow in a free electron gas, we shifted the Fermi 
surface in momentum space, thereby obtaining a super- 
fluid density Us = N and a phase stiffness '^ ep- If we do 
the same in a filled band we get nothing! A filled band is 
a single state with no degree of freedom, and shifting the 
distribution nk = 1 produces the same state. Hence the 
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obvious conclusion: a filled band has no stiffness. J can 
only result either from doping or from the order param- 
eter A that creates free carriers even if there were none 
to start from. As a result the phase fluctuation T^'^ is 
severely reduced: it becomes of order A, comparable to 
the BCS pair breaking Tc in every case. The difference 
with free carriers is dramatic, due to the narrow range of 
free carriers in momentum space. 

One more general comment is in order: is it possible 
to go smoothly from the semiconductor to the metal? Put 
another way, does a reduction of the gap Ao to drive us 
back to a regular Fermi liquid metal? As long as we were 
dealing with the density of states p only, the answer was 
"yes" in 2d: where free carriers sat in the Brillouin zone 
was of no concern. When calculating the phase stiffness 
J the answer is "no" . The standard situation is an indi- 
rect gap which closes when the minimum of the conduction 
band is lower than the top of the valence band. Thereafter 
free carriers exist even in the normal state, but they re- 
main localized near the corresponding band extrema. The 
insulator with no Fermi surface will not turn suddenly into 
a metal with a large Fermi surface (a sphere for free car- 
riers). Instead a small Fermi surface will appear near the 
band extrema, responsible for the very small superfluid 
density (and phase stiffness). Such a behaviour should be 
contrasted with the Id case, in which a smooth transition 
is possible. The standard example is the appearance of 
a commensurate spin density wave at zone edge: the re- 
sulting antiferromagnetic order produces Bragg scattering 
and opens a gap Aq. Near the gap the band dispersion is 
hyperbolic, returning to the normal metal behaviour when 
the distance (e — Eq) to the band bottom Eo exceeds Ao : 
in such a case the transition metal -^ semiconductor is 
smooth. But then the physics is completely different, as 
usual in Id geometries. We see no realistic 2d system that 
could provide a smooth transition and a large superfluid 
density (that would mean a gap spread over the whole 
Fermi surface and translating together with it). 

We now turn to the calculation of J, putting all the 
pairs in a state of total momentum Q and looking at the 
variation of the total energy. The BCS calculation for the 
energy is unmodified but for the following transformation 



C/c -^ ^k{Q) — (^k+Q/2 + C-k-|-Q/2)/2 . 



(30) 



where ^k is the dispersion relation of the band. This trans- 
formation changes only the kinetic energy: in our case we 
have for the two bands: 



KiQ) = 2J2vl[S^^{Q)-^,] + 2Y,vl[S^^{Q)-^,]. (31) 

ki^v kGc 



Since the shift of Q of a filled band does not change the 
kinetic energy of the band we can write: 



K{Q) = 2j2{a~t^) 

kG-u 



ken 



^k(Q)-M] 



kGc 



(32) 



We use at this point the hypothesis that in the region 
where Uk and Vk are significantly different from zero the 



two bands are parabolic. We can thus expand sz for small 



^k(Q)=^4--^ and S£(Q)^e^ + ^ (33) 



Sniv 



Srrir 



where rric/v are the effective masses of the two bands. The 
expression for K becomes: 

ken kec 

By introducing the 2-dimensional density of states per 
spin: p = m/2TT we obtain: 



K{Q) - KiO) = ^ 



-Aa P + OO 

DO JAo 



(35) 

and the masses simplify from the expression. We thus ob- 
tain for J: 



Ja 



TT 



2- 



e + /i 



e-A^ 



.M+pfT^ Vi^-i^r+^' 



(36) 

and since the integration is elementary we can write a 
closed expression for T^^ 



T, 



KT _ Vi^o - P? + A^ + ^{Ao + pr + ^' - 2/io 
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(37) 

We comment this result in the two cases of e = and of 
finite doping. 

In the undoped case the expression simplifies even more. 
In the superconducting region, Ao < Z\*, we find 



nKT 



(Zi„ - 2Ao)/8 



(38) 



T^ is thus comparable to the order parameter A, in con- 
trast to a BCS superconductor, a spectacular consequence 
of the severe reduction in phase stiffness. The behavior 
near the transition point is easily understood. There the 
gap is large as compared to A and the density of free car- 
riers vl lies in the tail of the BCS distribution: it is of the 
order pA'^ /Am- The stiffness is accordingly of the same or- 
der, leading to (|38|). Our conclusion is that the transition 
to the normal state is always due to unbinding of vortices, 
in sharp contrast to the weak coupling BCS limit. 

We consider now the case of large Ao- We can use the 
approximation used in (pfl) and keeping the lowest order 
in A/Ao in Eq. (|^ we find: 



rpKT _ 



1 



2AI 



(39) 



As expected T^^ is proportional to the stiffness of a single 
band with Fermi energy = e, this result is true for any sign 
oi p— Ao provided p ^ Ao 3> A. The factor of proportion- 
ality is greater than one indicating that the lower band is 
contributing to the stiffness. This result implies that even 
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in the doped system (for e — > 0) the lowest critical temper- 
ature is T^ due to binding-unbinding transitions of the 
vortices. The anomalous large value of the critical temper- 
ature found at mean field is thus never reached in reality. 
This result is not changed by a more realistic model for 
the dispersion relation as shown in Appendix ^. 



Conclusions 

Superconductivity is usually concerned with a regular metal, 
in which free carriers exist at the Fermi level. We show 
here that superconductivity can also occur in semicon- 
ductors as soon as carriers in either band feel an attrac- 
tion. The semiconductor is characterized by a gap 2Z\o, 
the origin of which is not specified. The attraction is char- 
acterized by the superconducting gap Am that would exist 
in a regular metal. The relevant parameters are the ratio 
Ao/Am and the doping of the semiconductor. 

Superconductivity may happen even in an undoped 
semiconductor in which no free carriers exist in the nor- 
mal state at T = 0, as soon as the gain in superconducting 
energy '^ Am exceeds the cost in producing carriers '^ Ag. 
Such a situation is completely equivalent to excitonic in- 
sulators, in which electrons repel instead of attracting: 
then electrons and holes attract and they form bound ex- 
citons. If the binding energy is larger than the gap these 
excitons appear spontaneously. In our model we find a 
sharp transition at a specific value Zi* ^ Am (the precise 
value depends on details of the density of states). The 
ground state is superconducting if Ao < A*, insulating if 
Ac > A*. Near the transition it departs from the familiar 
BCS regime: the order parameter A, the real quasiparticle 
gap Ag and the critical temperature Tc are such that 

Zi < Tc < Z\g 

That may be viewed as a "pseudogap" behaviour, mon- 
itored by the competition between semiconducting and 
superconducting features. 

In the undoped case only the vicinity of the transi- 
tion point Z\* is unusual. This is no longer true when 
the semiconductor is doped: then the sharp transition be- 
comes a crossover and superconductivity extends all the 
way. Nothing very spectacular happens when Ag < Z\* : 
the effect of a small doping is minor. In contrast the "semi- 
insulating region" Ao > A* is extremely unusual: that is 
the main result of our paper. Past the crossover the chem- 
ical potential first remains inside the band gap. We then 
find a superconducting state very reminiscent of the Bose 
Einstein condensation of preformed pairs. The order pa- 
rameter A is proportional to the square root of doping, 
and the characteristic energies display the same pseudo- 
gap ordering as in |2|. That behaviour prevails until the 
gap reaches A*'^ /e, where £ is a typical Fermi energy in the 
conduction band. Then the chemical potential enters the 
conduction band and, after a transient region, the system 
returns to a single band BCS superconductor that involves 
only extrinsic carriers. The range of anomalous behaviour 
is very broad, and thus generic. 



The main body of the paper is concerned with the sim- 
plest conceivable model, with a constant density of states p 
(a sensible choice for a 2d system), and a mean field treat- 
ment "a la BCS" , both at zero and at finite temperatures. 
The only new feature is the presence of the semiconduc- 
tor gap Ao- The problem is then studied both analytically 
and numerically. Both simplifications can be questioned: 
we now address the issues in succession. 

A constant density of states is not realistic for two 
reasons: (i) it ships the gap states to infinity, while stan- 
dard mechanisms for a small gap leave the states close 
to where they originated from, in an energy range ~ Aq 
(ii) in 3d the density of states p (^) starts as V? near 
band edge. Improvements are explored in Appendix H, 
where other densities of states are tried, still within a 
mean field picture. Numerical results show that the quali- 
tative evolution remains the same. In the undoped case the 
sharp superconductor-insulator transition is just shifted, 
the square root behaviour of A being unaffected. In the 
doped case the calculation can be carried out explicitly in 
the anomalous region, if 

A<^Ao-^i<.Ao 

Free carriers at T = are then mostly in the conduc- 
tion band and their density w^ is small, of order Z\^ . Our 
main result, A proportional to the square root of doping, 
is maintained (with modified coefficients). We conclude 
that the physics described in this paper is not too sensi- 
tive to details of p (f ) . We did not consider the case of a 
Fermi surface close to a Van Hove singularity, as found in 
a nested 2d system. That could definitely be done. 

Departures from mean field behaviour is a more deli- 
cate problem. Thermal fluctuations of the phase of A are 
controlled by the phase stiffness J. In ordinary metals J 
is comparable to the Fermi energy: phase fluctuations are 
small and a mean field picture is correct except very close 
to T*^^ (where fluctuations turn to a Kosterlitz-Thouless 
transition in 2d). In our semiconductor we show that J is 
only due to free carriers (it would be = for a flUed band) . 
J is comparable to the Fermi energy and thus small: it fol- 
lows that the phase fluctuation mechanism for Tc always 
prevails (at least in 2d and for Ao < A^/e): that is the 
second unexpected result of our analysis. Simple consider- 
ations are enough in order to provide orders of magnitude, 
but a real quantitative theory is still lacking. 

Granted these difficulties, we only tackled the prob- 
lem, our goal being to emphasize simple physical ideas. 
Many important questions remain open. For instance one 
should worry about the origin of the semiconducting gap 
Ao- If it arises from another instability (lattice distortion, 
magnetism, etc.), does superconductivity react on that 
primary instability? If it does we face a coupled problem, 
in which A depends on Ao and in reverse Ao depends on 
A. A self consistent calculation is required, and it may 
lead to unexpected flxed points ||I3]. Another issue is the 
relevance of interband pairing, as envisaged in |10 . Even 



if it is not the primary effect, it may affect the quantita- 
tive results. (It seems that such am interplay of intra and 
interband pairing is crucial in understanding the smooth 
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transition of the stiffness when the gap Ag closes to zero 
in a ID band.) These remarks open exciting possibilities, 
which we did not approach. As it stands our paper is only 
an exploration. 
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A Non constant density of states 

In this appendix we consider a more realistic form for the 
density of states, in which the states removed from the en- 
ergy region [— /\o, ^o] remain in the vicinity of the edge. 
We explore two simple cases: a discontinuous density of 
states (typical of 2d system) and a density of states that 
vanishes at Aq with a square root of ^ — Ag (typical of a 
3d case). These two possibilities are shown in Fig. |^. We 
solve the equations numerically showing that the qualita- 
tive behaviour found for constant density of states is not 
changed. 

Eqs. (ph and mu) can be easily generalized when the 
density of states is not constant: 



1.5 - 



1 - 



PiO 



v^(^)+fp[E{0] 



^-P 



tanh^^(0/2 '■+- 



d^ = 



PniiOd^ 






:d^ 



(40) 

In this equations PmiO and p(^) are the densities of states 
of the metallic and insulating phase respectively. The in- 
tegration path is now P' — [— w™,— Z\o] and [Z\o, -|-ijj„i], 
since we conserve the number of states by deforming p and 
not by shifting the whole band. The function p{^){> 0) 
is thus constrained by the conservation of the number of 
states when passing from the metallic to the insulating 
phase: 



PiOdC = 



PmiOd^- 



(41) 



It is clear that in order to fulfill Eq. (41) p{^) must depend 
on Ao too. 

Equations (Eoh are quite general. We will consider two 
specific simple cases. We assume that pm is constant, and 
that p(^) is defined in the two cases by the following ex- 
pressions: 



PliO/Prr 
P2{0/Prr 
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(42) 
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where x = {\^\-Ao)/Ao and p = for \^\ < Ao [cf. Fig. |I], 
the dotted line correspond to pi and the dashed line to p2\- 
The numerical results for A and /i in the zero temperature 
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Fig. 8. Zero temperature numerical results for the density of 
states p = pi (discontinuous one): (a) Chemical potential, and 
bottom of the upper band (dashed), (b) Order parameter A. 
ejAm take the following values: 0, 10~^, 10"^, 10~\ 1. 



case are reported in Fig. pi for the discontinuous density 
of states (/Oi), and in Fig. 19 for the second one. There is 
no qualitative difference with Fig. |[ Different shapes of 
the density of states correspond to different values of the 
ratio A* I Arm always keeping it of order 1^. 

Most of the results can be understood analytically. Let 
us first consider how the order parameter vanishes at the 
transition for no doping. In this case A is small and we 
can expand E in the second of (^). After the equation 
for the critical An is subtracted we obtain: 



= 



PAl{i,)-pAM) 



di 



^ Although the energy range over which p(^) varies is of 
order Ao, we have tried a more general form in which that 
energy scale is different from Ao. The maximal variation of 
A*o obtained is an increase of a factor e ~ 2.71 when all the 
removed states are accumulated at the edge. 
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Fig. 9. Numerical results for the density of states p 
(square root). Conventions are the same of Fig. H. 
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The last term for r ^ (Z\o ^ Z\*) yields simply A'^ir) 
multiplied by a constant. The second term vanishes like 
r^+" if p{^) ~ (C — ^o)"- This term is always subleading 
compared to the first term that vanishes linearly in r. 
We thus conclude that the dependence A ^ t^I"^ is quite 
general. This is not surprising for a mean field theory. 

For finite doping and T = it is interesting to study 
the dependence of A on e. In particular we can readily 
obtain this dependence analytically when Ao >• Am- We 
distinguish the two cases p < Ao and p, > Ao- When 
Ao ^ Ao — fi ^ A the lower band is nearly decoupled 
and the conservation of particles gives: 



SN = 2epr, 



P(0 



m-^^f 



A^ p{Ao 
2 Ao- 



m) 



/^ 



(44) 



We thus recover A ~ e^/^ for a general density of states, 
again a mean field exponent. Note that the coefficient in- 
volves the density of states on an energy scale ~ {Ao — p.) 



(measured from the bottom of the conduction band) . This 
fact does not invalidate the conclusion since the chemical 
potential is practically independent on e, determined as 
usual by the gap equation. 

For Ao large enough jj sits again in the valence band, 
and the limit Aq ^ p,— A^ A is reached. In this case the 
BCS theory applies in the upper band and one recovers 
the familiar exponential behavior of A on the inverse of 
the density of states at the Fermi energy. 

Concerning the critical temperature we have verified 
that using the density of states pi one recovers the same 
behaviour found for a constant density of states: near 
the transition point and for small doping it vanishes like 
1/w(1/t) or l/u'(l/e). In the BCS region {Ao » Al^/e) 
the usual BCS Tc is recovered. 

We can readily calculate the phase stiffness if we as- 
sume that the dispersion relation depends only on k^. In 
this case we can generalize the results for the variation of 
the kinetic energy: 

Sk(Q) =(,+a^+ ?^'^ -I- 0{Q^) , (45) 

where ^' = d^{k'^) / d{k'^) . Using this expression for the 
kinetic energy and the invariancc under rotation we obtain 
for each band: 



dE{Q) _ , ^ 

^ keband 
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(46) 



The derivatives of the dispersion relation can be related to 
the density of states: ^'f. — l/Airp and ^^' ~ -p'/Wir^p^. 
Putting everything together we obtain for our problem 
the same phase stiffness J as in (|36|), but with an ad- 
ditional factor A(^) in the integrand of (B6h [provided 

p{0^pi-0]-- 



A{0 



p'iO 



« 



p{x) dx . 



(47) 



When the density of states is constant A{^) is 1. One can 
show that the integral of [A{S) — 1] is equal to Ao'- it thus 
gives the contribution of the additional states to the phase 
stiffness. The various contributions to the phase stiffness 
are just rearranged and we find that the shape of p{^) 
does not affect qualitatively the dependence of the phase 
stiffness on the important parameters Aq and e. Moreover 
the quantitative change is not enough to let T^'^ become 
larger than the Tc given by particle excitations. Thus the 
phase fiuctuations should always be responsible for the 
transition to the normal phase. (Except for very large Ao 
at fixed £ where the BCS limit is recovered.) 



B Interband pairing 

In this paper we have considered pairing of time reversed 
state within the same band. If we indicate with aka and 
bka the destruction operators for the lower and upper 
band respectively, the possible pairings considered in the 
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main body of the paper are: (ak^a^ki) and {bi;]b^ki)- The 
two bands are then coupled by the interband interaction, 
which is assumed identical to the intraband interaction. 

Another model of pairing in a two band semiconduc- 
tor was proposed some years ago by Kohmoto and Takada 
]lOt (see also Q). The main difference with our proposal 
is that in that case the paired states are not time-reversed 
one of the other, specifically electrons of different bands 
are paired: {ak^h^ki)- (The spin was not included in the 
original formulation and it is introduced here for compar- 
ison with our work). The mean field calculations predicts 
at T = a first order transitions from an insulator to a 
superconductor. We discuss briefly in this Appendix this 
model [|0| and we compare it to ours in one simple case. 

The main difference arise from the fact that the kinetic 
energies of the paired states are different: one expects that 
the energy cost should accordingly be larger. Although it 
is possible to have interactions that favor only one of the 
two possible states, when the strength of the interaction is 
the same in the two channels the intraband pairing should 
be favoured due to the lower (kinetic) energy cost. We 
demonstrate this fact in a simple case. 

We consider a local interaction of the form 

-Ujd.pi.)p[.), 

where p{x) = Y.a^'^i^)'^'^ (^)o- ^^^d i'<y{x) = a^{x) + 
hij{x). The antiparallel spin component of this interaction 
can be reduced to two terms: 



Vi = -2U dx (a\al + b\b\] (a^a^ + b^b^) , (48) 

V2 = -2U fdx (a\b\ - a|6|] (b^a^ - b^a^) . (49) 

Vi would induce the pairing suggested in this paper, while 
V2 is needed in order to obtain the state proposed by 
Kohmoto and Takada. We consider the simplest case of 
constant density of states p, symmetric bands, and no dop- 
ing (^ = 0). We also keep the energy cut-off Ao + ojc = Wm 
finite. We can then evaluate the Hamiltonian on the re- 
spective trial wave functions. For the intraband pairing 
proposed in this paper the ground state wavefunction is: 

1^1) = n (^'^ + ^'^ ^I-T^-fci) ("fc + "fc "Ir^-fei) l«ac) 

while for the state proposed in Ref. ||10| the wavefunction 
takes the form: 



,«.n(i^)( 



1 - 4i^-fcT 



V2 



\vac) . (51) 



Here the singlet state has been chosen to take maximal 
advantage of the interaction. The energy difference with 
respect to the normal state (SE) is in the two cases: 



6E1 



2AoUJc + i^l - {Ao + LOc)^J{^o+^cy+^'^ 



+Ao^Al + A^ 

5E2 ~ pLOc [2Ao + LOc — 2UlUcP] 



(52) 
(53) 




Fig. 10. Comparison of the ground state energies of the state 
presented in this paper (full lines) and of the one proposed in 
Ref [|10| (dashed line). The upper part of the plot shows the 
order parameter for the two models. 



where the value of A in (p2|) is given by the usual gap 
equation: 1 — 2pU j." " di^j ^ i^ 4- Z\^. The minimum 
value of the interaction ([/""") for which 5E changes sign 
determines the transition from the insulating to the super- 
conductig phase. For the two states we have respectively 
pC/™" = l/21n(l -h ujclAo) and pf/a""™ = 1/2 + ^o/wc)- 
Thus C/™*" < C/™*" always. Moreover explicit evaluation 
of the energy shows that the energy gain of the intraband 
pairing state is always larger than the interband one. In 
Fig. |l^ we plotted bE in the two cases when bijJcjAo — 1. 
The transition is second order in case (1), first order in 
case (2). We report also the average of the order parame- 
ter over the states that participate to the pairing: 



Aav^ 



E.i 



(54) 



which clearly displays the nature of the transition. 



References 

1. H. Ding et al, Nature 382, 51 (1996). A.G. Loeser et al, 
Science 273, 325 (1996). 

2. P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 
195 (1985). 

3. See for instance: F. Gebhard The Mott metal-insulator 
transition. Models and method, Springer, Berlin, Germany 
(1997). A discussion on this point can be found in P. 
Nozieres, Eur. Phys. J. (in press). 

4. M. Randeria, Varenna lectures (1997) |cond-mat/9710224 
R. Haussmann Phys. Rev. B, 49, 12975 (1994). J. Ran- 
ninger, J.M. Robin, M. Eschrig, Phys. Rev. Lett. 74, 4027 
(1995). F. Pistolesi and G.C. Strinati, Phys. Rev. B 53, 
15168 (1996). S. Stintzing and W. Zwerger, Phys. Rev. B 
56, 9004 (1997). And references quoted in these papers. 

5. V.L. Berezinskii, Sov. Phys. JETP 32, 493 (1971). J.M. 
Kosterlitz and J. Thouless, J. Phys. C 5, 124 (1972). 



14 P. Nozieres and F. Pistolesi: From semiconductors to superconductors: a simple model for pseudogaps 

6. P. Monthoux and D. Pines, Phys. Rev. B 49 4261 (1994). 
J.R. Schrieffer, X.G. Wen, S.C. Zhang, Phys. Rev. B 39, 
11663 (1989). V.J. Emery and S.A. Kivelson, Nature 374, 
434 (1995). P.A. Lee and X.G. Wen, Phys. Rev. Lett. 78, 
4111, (1997). P.W. Anderson, J. Phys. Condens. Matter 8, 
10083 (1996). E. Dagotto, A. Nazarenko, A. Moreo, Phys. 
Rev. Lett. 74, 310 (1995). C. Castellani, C. Di Castro, M. 
Grilh, Z. Phys. B 103, 137 (1997). 

7. D. Jerome, T.M. Rice, W. Kohn, Phys. Rev. 158, 462, 
(1967). 

8. M.L. Cohen, p. 615 in Superconductivity, R.D. Parks Ed., 
M. Dekker, New York (1969). 

9. W. Hanke and M.J. Kelly Phys. Rev B 45 1203 (1980) 

10. M. Kohomoto and Y. Takada, J. Phys. Soc. Japan, 59, 
1541 (1990). 

11. J.M. Kosterlitz and J. Thouless, J. Phys. C 6, 1181 (1973). 
D. Nelson and J.M. Kosterlitz, Phys. Rev. Lett. 39, 1201 
(1977). 

12. V.M. Loktev, S.G. Sharapov, Con d. Matter Physics (Lviv) 



11, 131 (1997),|cond-mat/9706285| . E. Babaev and H. Klein 



ert, cond-mat/9804206 



13. A numerical analysis of the interplay between a supercon- 
ducting state and a Mott insulating state is given in: F.F. 
Assaad, M. Imada, D.J. Scalapino, Phys. Rev. B 56, 15001 
(1997). 

14. Y. Morita, Y. Hatsugai, K. Kohmoto, Phys. Rev. B, 53, 
8561 (1996). 



